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Abstract 

The coupled Bethe-Salpeter bound state equations for a QQ system, where Q = {U, D) 
is a degenerate, fourth generation, super-heavy quark doublet, are solved in several ladder 
approximation models. The exchanges of gluon, Higgs and Goldstone modes in the standard 
model are calculated in the ultra-heavy quark limit where weak 7, and contributions 
are negligible. A natural / = and 1 = 1 multiplet pattern is found, with large splittings 
occuring between the different weak iso-spin states when Mg, the quark masses, are larger 
than values in the range OATeV < Mq < O.STeV, depending on which model is used. 

Consideration of ultra-heavy quark lifetime constraints and U — D mass splitting 
constraints are reviewed to establish the plausibility of lifetime and mass degeneracy re- 
quirements assumed for this paper. 
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1. Introduction 

The prospect of SSC and LHC experiments to probe TeV phenomena has stimulated a 
wide variety of ideas about the physics at these energies. The standard model itself, expanded 
to a super-heavy fourth generation, could display some striking new effects. There is a variety 
of "heavy-neutrino" fourth generation models in the literature,^ and we adopt the attitude 
that some model of this type with the Higgs mechanism providing the fermion and gauge 
boson mass generation represents nature at 0.5-1.0 TeV. The effect that we explore in this 
paper is the strong Yukawa coupling of the Higgs sector to the fourth generation and the 
bound state spectrum in the heavy quark-anti-quark system (QQ) that is generated by 
this strong coupling. ^'^ There have been a number of heavy fermion physics studies since 
the Higgs mechanism and the standard model were proposed. ^'^ Like the heavy Higgs limit 
studies, the indications are that perturbative calculations break down and some new, strong 
coupling features of the standard model and/or totally new physics should emerge at the 
0.5-1.0 TeV energy ranges. 

In a previous work,^ we described some new, deep-binding features of the 0~ ground 
state of a single, super-heavy QQ system that were driven by the dominant, attractive Higgs 
scalar exchange plus weaker, but in some cases still significant, gluon exchange. The general 
feature was that a very strong binding compared to QCD alone was generally exhibited for 
Mq > 0.4 — 0.5 TeV in all of the model calculations that we tested. 

In this paper, we expand significantly on this theme by doing a coupled channel anal- 
ysis of a heavy quark doublet, {U,D), including the Goldstone boson exchanges as well as 
the Higgs exchanges, and we discuss the s-channel exchange as well as the i-channel ex- 
change contributions to the Bethe-Salpeter equation kernel. Here we include the results 
of calculations of the 0+ and 1~ masses in addition to the 0~ ground state case analyzed 
previously. 

Because of the weak iso-doublet structure of the super-heavy {U,D) system the bound 
states decomposes naturally into weak iso-triplet and weak iso-singlet states when Mu ~ M]j 
and Mz ^ M^. We find that the interplay between the scalar Higgs exchange and the 
(opposite sign) pseudoscalar Goldstone boson (longitudinal W and Z) exchange results in 
a clear iso-triplet/iso-singlet mass splitting in all of the channels in all model calculations 
we used. This signature is quite distinct for the fourth generation fermion doublet, since it 
must be aoDroximatelv degenerate according to oresent constraints from W/Z nhvsics. 



In the following section, we outline the coupled-channel Bethe-Salpeter formalism 
used to calculate bound state masses. In Sec. 3 we present results of the calculations in 
the different approximation schemes used and comment on the main features. In Sec. 4 we 
show that requiring meta-stability of fourth generation quarks puts only weak constraints 
on mixing with a lighter generation. We also review the evidence that a fourth generation 
of quarks, if they exist, would have to be approximately degenerate in mass. In Sec. 5 we 
summarize and conclude. In Appendix A, we describe the subtraction used to regulate the 
Bethc-Saleter equations. In Appendix B, we give a short analysis of the running of the 
Yukawa coupling and its relevance to the calculations presented. 

2. Strong Yukawa Binding of Ultra Heavy Fermion Weak Doublets: Bethe— 
Salpeter Equation Formalism Including Goldstone Boson Exchange 

If there are fourth and higher generations of fermion weak isospin doublets, the frac- 
tional mass-splitting within these ultra-heavy doublets^ is constrained by the standard, 
three generation model's remarkable survival after two decades of stringent experimental 
tests. A small mass-splitting, or near degeneracy, could lead to a classic isospin-multiplet 
pattern for the bound state spectrum, and we show in this section that this is indeed the case 
when the Yukawa couplings of the Higgs and Goldstone-bosons^ are incorporated into the 
bound state equations - ladder approximation, Bethe-Salpetcr equations in our study. We 
assume that the usual mass generation mechanism operates in the fourth generation sector as 
well, so the ultra heavy fermion doublets will have strong Yukawa couplings to the standard 
model Higgs and Goldstone modes. These Yukawa couplings, scalar for the Higgs boson and 
pseudoscalar for Goldstone bosons, become the dominant ones when the quark masses reach 
400-500 GeV (though the pseudoscalar exchanges actually become competitive only in the 
tight binding, highly relativistic regime. They decouple in the weak-binding limit). 

Wc develop below a set-up for the two, coupled-channel Bethe-Salpctcr equations of a 
quark doublet-antidoublet that interact through gluon and weak-sector Higgs and Goldstone 
boson exchange. 

Let us consider a doublet of heavy quarks, {U, D), interacting in the Feynman gauge 
by exchange of gluons, g, a Higgs boson, H, and Goldstone bosons x^)X°- The photon and 
and Z exchanges are always weak in the gauges chosen, have been checked and found to 
give insignificant contributions to the binding, and are not included in our development. It 
is convenient to adoot Fevnman eauee in the ladder aooroximation fnerturbative vertices. 



no crossed graphs), which we will use throughout this work. There are four quark-antiquark 
channels: IJU, tJD, DU and DD. These channels become coupled when the charged Gold- 
stone boson exchanges are included, as illustrated in Fig. 1, where the Bethe-Salpeter 
equations are shown graphically. As discussed below, these equations admit a subtraction 
at fixed q that eliminates the divergent, g-independent, annihilation graphs which enter the 
scalar and pseudoscalar bound state channels.^ The corresponding gluon annihilation graph 
does not require consideration because we study only color singlet bound states. 

The coupled-channel, doublet quark-antiquark bound state equations are shown pic- 
torially in Fig. 1 and displayed in Eqs. (2.1) and (2.2) below. The subtracted, decoupled 
versions of these equations, are presented later in this section, Eqs. (2.3a and 2.3b). Let 
us display, in shorthand notation, the structure of the ladder approximation, momentum - 
space Bethe-Salpeter equations for the bound state amplitudes x(gf, P) in Feynman gauge. 

U -U Channel: 

S-\q+)Xuui<l^P)S-\q-) = j -0^[-l,Xuuik,P)l.G^^''{k-q) 
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D — D Channel: interchange [/ <-> D in (2.1). 
U -D Channel: 

S-\q+)XuniQ,P)S-\q-)^l ^,[- l^Xuni^, Ph.G^" (k - q) 
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D — U Channel: interchange U '(-^ D in Eq. (2.2). 

In Eqs. (2.1) and (2.2), Xuuil^P) ^^'^■■> represent the Bethe-Salpeter amplitudes for 
the UU etc., systems with relative 4-momentum q and total 4-momentum P. The quark 
masses, m, are assumed to be deeenerate, q-l = a ± ^ and S-Ua+) are the inverse auark 



propagators with momenta q±. Here G^^{k — q) represents the gluon propagator, assumed 
to behave hke {k — q)~^ as | k — q |— > oo. The SU(2) couphng constant is denoted by g2. 

One may decouple the UU and DD equations by taking the combinations {tJU ± 
DD)/\^, and the resuhing Bethe-Salpeter equations, subtracted as described in Appendix 
A, for the bound state wave-functions Xud^Ii-^)-! X-{q.iP) ^^'^ X+{QjP) ^^ad as foUows 
iXoui^^ P) obeys the same equations as Xud)'- 
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(2.3a) 



'\m^ {k - qf - Ml 'aM^ {k - g)2 - M| 
The quark mass m = m,u = mjj is taken to be the same for U and D. Experimental support 
for this assumption is reviewed in Sec. 4. 

The important features of Eq. (2.3) are 

(i) X-1 XuD ^^'^ Xdu obey the same equation when Mz = Mw, and they form an isospin 
triplet of states. 

(ii) The Goldstone boson contributions to this x-, Xud ^ Xdu ^^ipl^t are, like the gluon and 
Higgs boson contributions all attractive in the ground state, 0~ system, (see Sec. 3 
below) 

(iii) the x+ state obeys a different equation whose Goldstone-boson contributions are all 
repulsive in the ground state, 0~ system. 

Eq. (2.3) and the fact that pseudoscalar exchange becomes comparable to scalar for relativis- 
tic systems combine to produce a clear splitting between the singlet and triplet of ground 
state bosons in the strong binding limit. The triplet is below the singlet in mass for the 0~ 
system. This is very reminiscent of a "tt ^ ry" situation, except that this appears to be purely 
dynamical with no spontaneous symmetry breaking, bound state Goldstone phenomenon at 
work. 

The iso-singlet, iso-triplet situation is reversed in the 0''" and 1~ spin-parity channels, 
as we describe in the next section. 



(2.36) 



3. Results of Relativistic Bound State Calculations of Iso-Singlet and Iso-Triplet 
Masses 

We employ two different approaches to solving the ladder approximation, Bethe- 
Salpeter equations for the = 0^ ground states of the weak iso-spin singlet and weak 
iso-spin triplet QQ system. We also report results in the same approximations for the 
= 1~ and 0^ states. The principle feature that emerges, as anticipated in the preceding 
section, is the splitting between the iso-singlet and iso-triplet states due to the pseudoscalar 
exchange (in Feynman gauge): repulsion in the iso-singlet case and attraction in the iso- 
triplet case in the ground state, 0~, channel. The iso-triplet pseudoscalar, ground state mass 
solutions are found to fall to zero for high enough quark masses in all of the calculations. 
The reverse is true in the 0+ and 1~ channels, and the isoscalar masses fall to zero at high 
enough quark masses in these states. 

We evaluate the bound state energies (masses) using both the covariant kernel, 

which we refer to as the covariant gauge ladder approximation, and the instantaneous ap- 
proximation to this kernel 

We have suppressed any reference to gauge-dependent spin structure in (3.1), and designate 
generic coupling constants and masses as g and M, respectively. The instantaneous approx- 
imation, so called because the coordinate space potentials are instantaneous, reduces the 
problem to three dimensions, yielding what are referred to as Salpeter's equations. These 
are coupled equations for positive and negative frequency amplitudes. One can solve the cou- 
pled system or set the negative frequency amplitudes to zero and solve the positive frequency 
problem by itself. 

Covariant Gauge Formalism 

In the covariant gauge ladder (CGL) formalism, we decompose the 0~ Bethe-Salpeter 
amplitudes into four independent spinor functions: 



Xp{P,q) =75 



(3.2) 



In our calculation, only Xo ^^'i Xi need to be retained, and the weak coupling relation 
y, = — y„/2m can be emploved to reduce the svstem of equations to a single equation for 



Xo- In Euclidean variables, we have 



{k - qy 

rn^gl 1 / o 3M| o\ /• d^k Xn , . 
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where M^=Mi^y=90GeV has been adopted since including weak boson mass difference af- 
fects answers only in the third decimal place, D = D{q'^ , M'^,cos6) = {—q'^ + — m^)^ 
-t-M^g-^ cos'^ ^, M'q — —P^ is the mass-squared of the bound state, 9 is the angle between 
P and q and — +1, —3 for iso-spin 7 = 1, respectively. As mentioned in the preceding 
section, the sign difference in the last term between the I — 1 and 7 = cases produces the 
extra attraction and repulsion compared to the purely attractive gluon and Higgs potentials. 
The method of solution is discussed in Ref. 3 in some detail. The resulting bound state 
mass values as a function of quark mass, m, arc shown in Fig. 2 in the curves labeled SI and 
Tl. In this calculation, the iso-singlet (SI) state's binding energy steadily increases, but the 
bound state mass never sinks to zero in the fermion mass range we have investigated. The 
more tightly bound triplet state (Tl) achieves zero mass when the quark mass reaches about 
900 GeV. Above 750 GeV quark mass, the sphtting between the different bound state masses 
becomes quite pronounced. The 0+ and 1^ calculations proceed in a similar fashion with 
the isosinglet being more tightly bound than the isotriplet. The details are not presented 
here, since they do not illuminate the discussion. 

Salpeter Equation Formalism 

Making the approximation 3.1b, integrating over qo in the Bethe-Salpeter equation, 
and projecting out the 0~ channel amplitude, one obtains coupled equations for the positive 
and negative frequency amplitudes, X±(?), where g =| ^|. The equations read 



- 2u;)x+(9) = ^ / dq'q'(v+x+ + V-X-) 
(e + 2u;)x_{q) = — / dq'q'(V-X+ + V+x.) 



(3.4) 



where lo — \Jq^ -\- rm?, E is the bound state eigenvalue and 
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The coefficients Cy and variables Zy etc. are summarized in Table 1. Details on the method 
of solution are given in Ref. 3. 

The 0~, ground state masses as a function of quark mass is shown for the positive 
frequency only case in Fig. 2, curves S2 and T2. Again the splitting between isosinglet (S2) 
and isotriplet (T2) bound state masses is striking for the heavy quark mass region, and the 
iso-triplct system becomes ultrarelativistic in the region above m — TSOCeF, producing a 
zero mass bound state at m ~ llOOGeV^. 

The solution for the positive plus negative, fully coupled, system is much more tightly 
bound as is seen in Fig. 2, curves S3 and T3, where a dramatic departure between iso- 
triplet (T3) and iso-singlet (S3) masses sets in already at m = AQ{)GeV , and the iso-triplet 
mass plunges to zero at m = 5206*6^. Figure 2 also shows the feature of the instantaneous 
approximation that differs from the CGL approximation, namely the turn-over (and eventual 
fall to zero) of the iso-singlet bound state mass as a function of the quark mass. This effect 
does not yet appear in Fig. 2 for the iso-singlet positive-frequency-only case, but it does 
occur at a quark mass above 1.20 TeV for this case as well. 

Focusing on the iso-triplet ground state mass values vs. quark mass curves in Fig. 2 
for the three solutions which were discussed above, we see that the general features agree 
though the different relativistic bound-state approximations produce different bound state 

mass values for a given quark mass. One can conclude that the Goldstone boson Higgs 
exchange plays an important role in the calculations, that deep binding at or above 500 GeV 
quark mass occurs, and that a dramatic weak iso-spin mass splitting is produced in the QQ 
spectrum of a degenerate, or nearly so, [/, D doublet system. 

With the role of iso-singlet and iso-triplet reversed, the features just outhned are 
present also in the = 1~ and 0''" bound states, and we display results of calculations for 
this svstem in Fies. 3 and 4. The nrinciDle distinction between 0^ on one hand and 1^ 



and O"*" on the other is the iso-singlet/triplet sphtting reversal and that the latter two are 
typically less tightly bound. 

In the next section we take up the question of ultra-heavy quark decay lifetime and 
intra-doublet mass splitting constraints. 

4. Constraints on Super— Heavy U, D Quark Lifetimes and onU — D Splitting 

Two related questions arise when one considers the phenomenology of super-heavy 
fermions. A crucial consideration for bound state physics is the comparison of quark lifetime 
to the period of bound state motion. The situation for the heavy top quark case has been 
considered by a number of authors/^ and we apply the argument of Strassler and Peskin 
to the super-heavy, fourth generation in this section. Related to the lifetime question is 
the U — D mass-splitting question. If the splitting is larger than the H^-mass, then the 
higher mass quark can decay to the lower mass quark by direct VF-emission, presumably 
with no Cabbibo, Kobayashi-Maskawa (C-K-M) suppression. The width would then be too 
broad to permit formation of narrow bound states containing the heavier quark. We review 
constraints imposed by measured W and Z-boson properties on the U — D mass splitting, 
which we would like to take as degenerate to a first approximation. 

Lifetime for decay to light quarks 

The lighter of the C/, D doublet partners must decay to lighter generation quarks, q, 
and we assume that one transition is dominant. Calling the corresponding C-K-M factor 
VQq, we have 



Mr 



3 



rQ^,~(180)Mel-[^J \Vq,\' , (4.1) 

where Mq » Mq^M^ is assumed. To make a conservative estimate, we will include 
only the gluon binding in estimating the time needed for bound state formation in a non- 
relativistic, weak binding approximation. Including the strong binding, Yukawa interactions 
will only improve the prospects for bound state formation. For pure QCD coupling, the 
characteristic radius is oq = (|q;sMq/2)~^ where is the effective QCD couphng strength 
for the problem (with a characteristic velocity Vn ~ {A/3)as/n for the n*'* radial excitation). 
Taking the ratio of twice the diameter to the velocity as an s-state formation time,"^^ one 
has t/orm ~ 2a'H4 ^ ~ 1, 2 ■ ■ ■ and in the ground state t/o^-jn ~ 9/2q^Mq. The inverse 
Bohr radius is the appropriate scale at which to evaluate ctg for the problem, and using the 
one-loop parameterization for ag 



where d — 12/(33 — 2nf), along with ^ = ^|as(/x^)MQ^ , leads to a condition on q;s(/x^): 

We show a plot of Q;s(l/ao) as a function of Mq/Aqcd in Fig. 5. As expected, since the 
binding scale is much less than the value of Mq, "^(^^^ is significantly larger than Q;(Mg). 
The latter is also displayed for comparison purposes in Fig. 5. 

Equating the time of formation and the decay lifetime, one obtains an upper bound 
on the value of | Vq^ | for given values of h-QCD and Mq. Choosing A-qcd — 0.20Ge\^, and 
Aqcd = 0.10 we show | Vqq \max vs. Mq in Fig. 6. 

Fig. 6 clearly shows the point, a surprise to us, that the constraints on Vq^ are 
rather weak, permitting for example Vqq ~ 0.3 ~ sin^c — Vus for Mq ~ b^QGeV with 
^QCD = 0.20G'ey. As mentioned above, including the Higgs boson, Yukawa coupling effects 
will only improve chances of bound state formation and further weaken the constraints on 
{VQqjmax- This observation is encouraging for prospects of detecting narrow, deeply bound 
states at higher energies. 

Mass splitting between U and D 

Turning to the question of mass splitting between the U and D quarks, we use the 
definition^^ 

and the parameterization of radiative corrections in terms of the factor Ar to write^^ 

Ml f,_Ml)^ ™ 1 1 (4.4) 



M| V M| y ^/2GF M| 1 - Ar 
The quantity Ar can be expressed as 

Ar = 0.071 - '-^^SpH , (4.5) 
sm Ow 

where the heavy fermion factor Spu is given by 
and 

AM^ = M? + M|- tf^S/ nHi . (4.7) 

^ ^ M? - M^ Mo ^ ^ 



The factor Ci — 3 for super heavy quarks and Ci — 1 for super-heavy leptons.^^ The masses 
are designated as Mt for the top quark and Mi and M2 for super-heavy fermion generation 
members. Using the 1992 data book^^ values for M\y and Mz and the value sin^ 9w — 0.230, 
we find from Eqs. 4.4-4.7 that the one standard deviation constraint is 
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For the case of one super heavy-quark generation and a degenerate super-heavy lepton 
generation, the plot of the bound on | Mxj—Mjj \ vs. Mt is shown in Fig. 7. There is evidently 
room for considerable mass-splitting in the bound (4.8), where | Mu — Mjj |< lAQGeV for 
Mt — ISOGey , for example. The dependence of the bound on the average value of the heavy 
quark mass is very weak and can be ignored for our purposes. 

An independent bound can be obtained from Z decay data. The leptonic and hadronic 
decay widths, for example, provide useful constraints. Defining the quantity sin^ 9w = 
k{Mz) sin^ 9w, one has approximately^^ 



sin 9w 
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with k{Mz) = 1 + (^3Q;/167rsin^%^ (^AM|/M|^. The value of a{Mz) is obtained from 
electroweak corrections to a that do not include the heavy quark effects. Using the value 
Ar = 0.071 to define a{Mz) and adopting the first factor in brackets in Eq. (4.9) as the 
value of sin^ 9w to be used in the expression in the second bracket, we find 



sin^ 9w = 0.235(1 - 3.50 x 10~^^^) ■ 



(4.10) 



The expressions for the decay widths are 
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The factors p and 5ft in Eqs. 4.11 are given by 
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(Mf /M|), where sin^ 0w 



«(! + 4 sin^ %/3 - 16 sin^ 0w /9) /Stt sin^ 0^ cos^ 0^ 
0.235 is used from Eq. (4.9) onward to evaluate all of the above expressions. The constraints 

on I Mjj — Mjj I provided by the 1992 Particle Data Group vahies of r{Z and 
V{Z ^ hadrons) are illustrated in Fig. 7. We see again that significant mass splitting in 
a hypothetical fourth generation of fermions is allowed. The constraints on splittings are 
essentially independent of average mass. 

To conclude this section, we note that our assumption of degeneracy between super- 
heavy generation fermions is consistent with the present bounds, but that there could be 
shifts in our iso-triplet vs. iso-singlet QQ bound state masses due to ultra-heavy quark 
non-degeneracy effects. The pattern of a mass-splitting between iso-singlet and iso-triplet 
states of a given spin and parity would remain, however. 

5. Discussion and Conclusion 

In the highly relativistic, strong binding regime, we have shown that the contributions 
of the pseudoscalar, Goldstone-boson degrees of freedom to the binding of ultra-heavy QQ 
systems play a crucial role in determining the bound state spectrum. Even though these 
pseudoscalar interactions become negligible compared to the gluon and Higgs exchanges in 
the weak binding limit, they are solely responsible for the large iso-singlet vs iso-triplet 
QQ bound state splitting which we found for ultra-heavy, fourth generation, Q = {U, D) 
quarks. Depending upon the specific approximation scheme used, we saw (Figs. 2,3,4) 
that the splitting between / = and 1 = 1 QQ states become large for Mq in the range 
QATeV <Mq< Q.STeV. 

Whether a splitting between bound states is large depends of course on the decay 
width of the (quasi) bound states, and our estimates of | Vg^ | in Sec. 4 and the analyses of 
decay widths of similar states into qq, H, Z etc. in Ref. 17-21 indicate that the splittings 

between states will be large compared to the decay widths for the Mq values in the range 
given above. The states should, therefore, be clearly separated. 

For the production and decay of neutral QQ (and lepton-antilepton) ultra heavy bound 
states, Refs. 2 and 17-21 provide an encouraging picture of the production and decay signals 
for such states. The binding indicated by our relativistic, bound state calculations is much 
stronger than that considered in Refs. 17-18 (gluon non-relativistic potential model only) or 
3. 19-21 fHiees non-relativistic model or Hiees dIus eluon non-relativistic Dotential model). 



so the prospects for detecting new heavy bound states, should a fourth generation exist, are 
made even brighter by our results. 

Avenues for further work include the detailed phenomenonlogy of the iso-multiplct 
system, in particular the possible production and decay of the charged state, and the ap- 
plication of these techniques to TcV scale supersymmetry bound state physics. We plan to 
address these issues in the future. 
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APPENDIX A 

We describe here the substraction procedure used in regulating the B-S equation. 
Referring to Eqs. (2.1) and (2.2) in the text, we see that the direct, P^-pole terms in each 
equation are independent of q. All of the other terms in the integrands depend on A; — g and 

they decrease like [k — q)~'^ ~ q~'^ as | g |^ oo for fixed k. Therefore, if the integrals exists, 
one finds that for fixed P, taking (2.1) in particular, 

. 92 i rr. ( [ d'^k ^ 

^ 4 M^P^-Mjj [J (27r)4^^^^'''^^ 

(Al) 

If S'~^(g_(_)|g[^oo ^ ^ as is the case for the perturbative propagators which we use in this 
paper, then the fixed P, large q limit gives 

xiq,P)^^ . 

In this case, the integral / d'^kx{k, P) does not exist and the integral equation has no solu- 
tions, in accordance with direct numerical studies which we have carried out. 



If the g-independent direct (p)-channel pole terms are absent, however, then the 
dependence of the integrands of all of the other terms in the large limit means that 
x{q, i-*)|g|^oo ~^ behavior is consistent with the existence of solutions, and we find that 
this is indeed the case. 

To render finite the equations (2.1), (2.2) and their counterparts with U and D in- 
terchanged, it is sufficient to subtract each equation at a fixed say g = 0, because this 
subtraction removes the ^'-independent, divergent term. To make this subtraction system- 
atic, one regulates {i.e. cuts-off) the integrals, makes the subtraction, and then removes the 
regularization. The subtracted equation, is as follows: 



S-\q+)Xuu{<l^ P)S-\q-) - S^' (^)xuui^, P)S-' y ^ 
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(A.2) 

If the first term on the left hand side of (A. 2) is equal to the first integral on the right hand 
side for each q (i.e., a solution for each q is found), then the companion terms, where ? = 0, 
on each side are guaranteed to be equal and the result is an eigenvalue equation for P and 
q. This equation is, of course, just the original equation without the offending, divergent 
term and we use this form as our regulated, finite Bethe-Salpeter equation which is to be 
projected into states of definite quantum numbers. 



APPENDIX B 



Running of the Yukawa Coupling — The Landau Pole 

In this appendix we briefly discuss the effects that arise because of the momentum 
dependence, or running, of the Yukawa couplings. We work at the one-loop level to see 
where the Landau pole occurs when the quark masses are in the deep binding region that 
we probe. O.STel/ ^ Mn ^ l.OTeV. 



The one-loop renormalization group equations for the gluon and Yukawa couphngs, gs 
and qy, without the weak gauge interaction effects, read^^ 

d 1 



,3 



and 

d 1 /9 



[B.la) 



where d = 12/(33 — 4np), Ug is the number of generations, and t = with /j, 

arbitrary renormahzation scale. 

Parameterizing the solution to (B.la) in the standard fashion, 

ATT^d 27r^d 



an 



2/A2 \ t 



where Aqcd — 0.2 GeV and /i — Aqcd is adopted, we can rewrite (A. lb) as 

dgy _ 9 o d 

-dF - 32^^- -1^- ■ ^^-^^ 

The solution to (B.2) can be written in terms of an integration constant C, to be fixed by 
the boundary condition: 

16^ l-2d 

Checking several limiting cases of (B.3), we see that 
when a^pr, = o? = 0, as it should. The characteristic Landau pole appears at t = ^ — 



In the weak Yukawa coupling limit, which is the light quark case, where only the Qg effect is 
kept, 

9Y{t) = , 

so the quark mass, proportional to gy and the weak scale v, behaves as 

V 1 

which is the well known leadine-loe, asvmntotic behavior of the lieht auark masses in QCD.^^ 



The constant C in Eq. (B.3) can be fixed in terms of the input quark mass, Mq, by 
requiring the condition^^ 



,2 



Ml^gl{2MQf- , (BA) 



for example. The requirement (B.4) can be solved for C to give 



with 

2Mq 



tq = in- 

J^QCD 

The expression B.3 with C evaluated by (B.5) can now be used to assess the large \ q\> Mq, 
behavior of gyit)- We can use the results (B.5) and (B.3) to determine Cq and qpoie, the 
location of the singularity for various Mq values. There is no longer a pole at large t, but the 
coupling still blows up at large t. Interesting values are: Mq = 0.5TeV, where Cq = —0.282 
and qpoie = ISTeV; and Mq = l.OTeV, where Cq = —0.288 and g^o^g = 3.5TeV. The t = 
pole in (B.3) is an artifact of the standard parameterization of qg, which is only applicable 
when I g |>> Aqqd- Qy should be sensibly constant for small momentum transfers. 

The momentum values that are important in the wave function in the Bethe-Salpeter 
integretion are of the order of the inverse Bohr radius, as a rough rule of thumb. Only in 
those cases when Mq k> ITeV and in the hmit as the bound state mass plunges to zero 
does the gy singularity come near to the relevant range of integration in the Bethe-Salpeter 
equation. Consequently, the binding energies as a function of Mq should not be substantially 
affected by the gy singularity below the point where the bound state masses fall to zero, 
even for the Mq ^ ITeV cases. 



Figure Captions 

Fig. 1: Bethe-Salpeter equations for the UD and UU bound states. 

Fig. 2: Isosinglet and isotriplet masses for mesons in the 0~ channel. Curves SI (isosinglet) 
and Tl (isotriplet) are the results from covariant gauge formalism. S2 (isosinglet) and T2 
(isotriplet) are the results from Salpeter's equation using only positive frequency components. 
S3 (isosinglet) and T3 (isotriplet) are the results from Salpeter's equation using both positive 
and negative frequency components. 

Fig. 3: Isosinglet and isotriplet masses for mesons in the 1~ channel. Solid curve (isosinglet) 
and dotted curve (isotriplet) are the results from covariant gauge formalism. Short dash 
curve (isosinglet) and long dash curve (isotriplet) are the results from Salpeter's equation 
using only positive frequency components. Short dash-dot curve (isosinglet) and long dash- 
dot (isotriplet) are the results from Salpeter's equation using both positive and negative 
frequency components. 

Fig. 4: Isosinglet and isotriplet masses for mesons in the 0+ channel. Sohd curve (isosinglet) 
and dotted curve (isotriplet) are the results from covariant gauge formalism. Short dash curve 
(isosinglet) and long dash curve (isotriplet) are the results from Salpeter's equation using 
only positive frequency components. Dash-dot curve is the isosinglet result from Salpeter's 
equation with positive and negative frequency components. The isotriplet result for this 
case is almost identical to the corresponding result obtained by including only the positive 
frequency components in the Salpeter's equation. 

Fig. 5: The value of as (solid curve) relevant for the bound state Bethe-Salpeter equation, 
ignoring higgs interaction, as a function of the quark mass. The dashed curve shows q;s(Mq) 
which is significantly smaller than the value of ag that should be used for bound state 
calculations. 

Fig. 6: The maximum value of the CKM matrix element | Vq^I allowed in order for the fourth 
generation quarks Q to form a meson bound state, q being a lower generation quark. This 
constraint is derived by including only QCD effects. Inclusion of the Yukawa coupling will 
allow even higher values of \VQq\. 

Fig. 7: The maximum allowed value for the fourth generation quark mass difference \Mu — 
M[)\ as a function of the top quark mass M^. The three curves represent the constraints 
using p parameter (solid curve), T{Z — > hadrons) (short dashed curve) and r{Z — > Z+/^) 
(long dashed curve). 
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Table 1 

Coefficients for the 0~ channel Salpeter equation kernels, Eq. (3.4) and below. Positive signs 
indicate attractive interactions and the negative sign a repulsive inter ation. 



